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ABSTRACT 



This investigation was undertaken to determine the feasibility of 
computing the friction velocity from available micrometeorological data. 
This could be of great assistance in finding the wind velocity profile by 
Deacon’s equation and the eddy shearing stress in the lower portion of the 
friction layer. 

This investigation was carried out at the U. S. Naval Postgraduate 
School, Monterey, California, during the period February 1959 - May 1959, 
in partial fulfillment of the requirements for the degree of Master of 
Science in Meteorology. 
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1. Introduction 



The friction velocity, a stability parameter, and a term related 
to the Richardson number word computed using micrometeorological data. 

The equations used for the computation of these three parameters are 
derived. These equations are unique in that they combine the different 
theories of two well-known Investigators of surface layer turbulence into 
an equation for the friction velocity and an equation for the wind profile 
in the surface layer. 

The micrometeorological data used in this Investigation is discussed, 
as well as the instruments used in obtaining the data and the observa- 
tional errors. A theory is presented on errors of aerated temperatures 
close to the ground. 

After modifying the derived equations into a form suitable for using 
with an electronic digital computor, the so-called "constants" in these 
equations are discussed. The advantages and disadvantages of using an 
electronic computor for this type of an investigation are then summarized. 

Various figures and graphs are used to Illustrate the diurnal varia- 
tion as well as the possible variation with height as determined from 
this investigation. A comparison of these results with those of previous 
Investigators is made, indicating close agreement in most eases. 
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2. Derivations 



Deacon [2], in the derivation of hia wind profile equation, as well 
as most other well-known investigators in the field of surface layer 
turbulence, considered only two virtual turbulence parameters - namely, the 
mixing length ^ and the friction velocity u A , where 




Here is the eddy stress and p is the air density. Lettau [10] introduces 

★ 

a third parameter, the vertical mixing velocity w, that is quite Important 



if nonadlabatlc rather than neutral conditions are Investigated. He also 



has a term which he calls the friction velocity u, where 



** o- 

uw * _x 

/° 

However, when the friction velocity is mentioned in this paper it will refer 
to the first definition. 

The surface layer is defined as that layer, immediately above the 
earth's surface, in which the eddy stress is constant with elevation through- 
out the layer. Eqs. (1-9) are descriptive of this layer. The notation^ is 
that employed by Lettau in [10]. 

In a neutral surface layer the equations below are applicable. 



★ * 
u * w 
a a 



( 1 ) 



^Throughout this paper u will be understood to mean the wind velocity 
at some level averaged over some Increment of time - say six minutes. 
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( 2 ) 



Z 



a 



k(z + z ) 
o 



k * von {Carman's constant 
z ^ m roughness parameter 



★ 

U 

a 




(3) 



★ ★ 
u v 
a a 



2 * 



(4) 



a 



* 



a 




(5) 



It should be noted that In Eqs. (1-5), u , w , and are constant 



with elevation in the neutral surface layer. 



In a non-neutral surface layer the next four equations are used. 



★ ★ 
U i w 


(6) 


** 

uw * _£ 

P 


(7) 


u * -if 


(8) 


T 


(9) 



In [10], Lettau derives a relationship between w and w^ which is of 
fundamental importance in the development of the remainder of this paper. 
Two basic assumptions are necessary in the derivation of Lettau 's relation- 
ship: 

it 

1. The ratio w does not depend upon variations of thermal stratifi 
cation, so that 



★ 

v 



~F 



* 

wa 



T~a 



( 10 ) 



3 



I 



Eq. (10) indicates that the vertical mixing length is changed in exactly 
the same ratio as the vertical mixing velocity if the thermal stratifi- 
cation varies. 

2. The second assumption is that the buoyant acceleration influences 
only the vertical parameters of turbulence. This assumption is necessary 
in order to fix absolutely the variation of mixing velocity and mixing 
length with thermal stratification. In an adiabatic atmosphere the term 

can be considered to represent the vertical component of a turbulent 
acceleration. Thus we have the vertical acceleration equation: 



ture with respect to height, and * mean absolute temperature averaged 
with respect to height In the layer considered. Eq. (11) Indicates that 
the mean vertical turbulent acceleration In a non-adlabatlc surface layer 
equals the adiabatic acceleration plus the buoyant acceleration. 




(ID 



where g * acceleration due to gravity, b0/dZ * change In potential tempera- 



Dpon multiplying Eq . (11) by 2 





but from Eq. 




, so upon substituting 




multiplying through by 




gives 



4 



or, 




tM 

T m &* 




rj§4 2 



Let 



x = 






x 






( 12 ) 



Therefore, x = — / 

j£ 



(13) 



Thus we finish with Lettau’s derivation by putting 




in terms of x. 



A = 



i_ 

(i+x) 



(14) 



The parameter x is important in the computations to be made in this paper. 

The next derivation was originally developed by Martin [12] combining 
Deacon’s wind profile with Lettau’s vertical mixing velocity. 

Starting with Deacon’s Equation 



in = U* / * + (15) 

iZ R2 a ^ 2„ / 



where /3 , as defined by Deacon, is a decreasing function of the Richardson 
number. Combining Eq. (7) and (8), one obtains 
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★ 

v^bere u 



X 



CKA 

az 



* 

and w 




So that 



u . = 




dm 

az 



Upon substituting into Eq . (15) 





an 

~5z 







But from E 4 . (14) ^ = I 

A ( 1 +x ) 

the radical by A/A 



, so that after multiplying inside 



_ 

bZ 




A 

l+X) o+*7 



/? Z* 





Dividing through by 




i/IK - & &L /z±zs.)~^ 
V az lez^o+x) V 2 <> / 
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Squaring out: 



lid = /ztz o r 2/3 

az pzfa+x? (. z, / 



But from Eq. (2) — I? ( +■ Z 0 J , so that 



<^U fe Wl (Z.'t'Zo) / 2 +g 0 1 ^ 

dZ k x z2(j + X)*- \ / 



or, 

= ^a. /ztZo l /zt2. ) ~ 

<)Z kz o 0+x) z l Zo/( 



Further combining of terms gives the equation in its final form. Namely, 



C>u _ A- [z^±Z° )'~ 2/3 < 16 > 

az iez ( ,(i+>') z ( z <> / 



Thus Martin has eliminated u + from Deacon's Equation - a parameter that 
is difficult to determine with accuracy. 

Rewriting Eq . (16) in the form 



<^u — 


fZ J r1o\ 


1-/3 


I /zf-2»r^ 


az 


Ty^W v ) 




iez 0 ( z* / 



( 17 ) 
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and comparing it with Eq . (15), one can easily see that an equation for 
the friction velocity has evolved: 



Wy / 2 + Z„\ ' 
(/+ y) 2 \ z o / 



(18) 



Using this equation to find one needs only to know the adiabatic mixing 
velocity, roughness parameter, vertical temperature gradient, and the 
stability parameter. 

This completes the necessary derivations. The remainder of the paper 
deals with the solution of Eq . ( 2) for x, Eq . (16) for ^3 , and Eq » (18) 
for u^, using available micrometeorologica 1 data. 
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3. Data 



Immediately upon deciding to conduct this investigation into micro- 
meteorology a search was initiated to find suitable data. It soon became 
apparent that there was a scarcity of micrometeorological data of the 
caliber necessary to conduct individual research into surface layer turbu- 
lence. Several sets of available published data had to be discarded 
because one or another necessary parameter had not been observed, the period 
of observation was too short, or the heights of the observations were not 
at suitable locations. 

The best data that could be found within a reasonable time was that by 
Gerhardt [5]>in 1949 and by Gerhardt [6] in 1950. The 1949 data that was 
of interest in this investigation consisted of wind observations at 12* , 

41*, and 91 1 ; aerated temperatures at 3', 6', 12', 20 1 , 35', 55 1 , 80', and 
110 1 ; non-aerated temperatures at 1 M , 6", 18", 3', 5’, 9', and 15'. The 
1950 data included in addition to the above, wind at 3 1 and non-aerated 
temperatures at 23', 33’, 45', 60', 7 8 1 , and 99’. 

The information was gathered during runs of from 24 to 44 hours for 
different seasons of the year, with an average run being about 30 hours. 
Observations were taken every half hour during a run, with the final value 
for that time and height being an average of five readings taken over a 
six-minute interval . 

Jehn [9] and Gerhardt [7] give a detailed explanation of the instruments 
used and their accuracy, the topography of the micrometeorological site 
area, and the vegetation at the site and surrounding fields. The wind 
measurements were taken with Bendix-Friez Aerovanes to the nearest mph. 

The errors involved were less than 1 mph for wind speeds in the range 
10-50 mph, increasing to nearly 2 mph after 13 months of operation. For 
wind speeds less than 10 mph the errors were slightly greater, and readings 
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below 3 mph were unsatisfactory. 

The wind measurements, because of their inaccuracy and the spacing 

*v 

of instruments, were a severe handicap in this investigation. However, 
it was hoped that by analyzing a large number of observations, the errors 
would be somewhat minimized. 

The temperatures were obtained from calibrated rod-shaped thermistors 

(ceramic resistors), about 0.75 inches long and 0.03 inches in diameter. 

All thermistors were fitted to a single calibration curve giving an 

' O o + o 

accuracy over the range of 0 C to 50 C of - 0.1 C. The aeratec} thermistors 
were mounted in anti-radiation housings, being ventilated by a conven- 
tional motor and exhaust fan arrangement at the top of the housing. The 
» 

non-aerated thermistors were coated with an anti-radiation paint and, in 
later runs, were shielded from both solar and terrestrial radiation without 
reducing the natural ventilation available. 

Because of all the precautions taken in calibration and shielding, 
one could expect that there would be little or no difference between the 
aerated and non-aerated temperatures. However, during this investigation 
it soon became evident that there were large differences (in some cases 
nearly 1°C) between the two elements, and in several instances one gave an 
increase in temperature with height while the other gave a decrease over 
the same layer* In order to determine which would be best suited for this 
investigation, a thorough analysis of both sets of temperatures was con- 
ducted for two runs. It was found that at elevations greater than five 

feet the non-aerated temperatures were higher than the aerated both day and 

o 

night, with the differences usually .2-. 6 C,for all pairs of thermistors. 

For elevations lower than five feet, however, it was discovered that the 
non-aerated thermistors were cooler during the early afternoon and warmer 

during the early morning, with the differences increasing with decreasing 
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elevation. This at first seemed puzzling, since it is contrary to what 
one would intuitively expect -- at least for early afternoon. However, 
the author believes that this can be reasonably explained. 

The magnitude of the vertical temperature gradient decreases rapidly 
with height near the ground, especially during early morning and afternoon. 
Now consider the ventilation system drawing in an equal volume of air from 
both above and below the level z. Because of the difference in tempera- 
ture gradient the absolute value of the temperature difference between 
z and z -az will be greater than that between z and z + ^z. This dif- 
ference will be greatest during early morning and early afternoon, and 
thus the aerated temperature will be warmer than the true temperature at 
z during the afternoon and cooler in the morning. 

Because of this effect, it vas decided that the non-aerated thermistors 
would give a more representative vertical gradient. However, of the four 
runs used in this paper (Runs 8, 10, 16, and 17) only on Run 17 was it 
possible to use the non-aerated temperatures since on Runs 8 and 10 they 
were only observed to 15* and on Run 16 the 3 1 element was broken. Even 
so, it is felt that greater accuracy in the temperature measurements are 
not warranted because of the large errors involved with the wind 
velocities . 

Gerhardt, in his final evaluation of the Micrometeorological Research 
Project (from which the data used in this investigation was obtained), 
stated: 

The importance of an accurate wind gradient, particularly in the 
lowest layers, cannot be overestimated; 12 measurement levels below 
50 feet would not be excessive. Great care should also be taken t q 
obtain instruments with as nearly identical response as possible. 
Although thermistors have certain undesirable characteristics, it is 
believed that the installation as set up was satisfactory, and that 
accurate and representative air temperatures were obtained. 
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4. Computations 

In this section, the working equations using finite differences in 
place of derivatives are illustrated. The evaluation of the supposedly 
"non -variable'* parameters is discussed, and the advantages and dis- 
advantages in using an electronic digital computer for this type of com- 
putation is pointed out. 

The particular working equation used to compute x was obtained by sub- 
stituting Eq. (2) into Eq . (12) and substituting finite differences in 
place of the partial derivative . This results in the working 

equation 



X = 






(19) 



The terras inside the bracket are constant for the layer during most 

★ 

of the run, although a new value of w was selected every 12 hours. It 

is easily seen from this equation that x is zero for neutral conditions, 

greater than zero for unstable conditions, and less than zero for stable 

2 

stratification. Values of (1+x) were tabulated. 

The next variable to find was ^3 . The partial derivative was 

first calculated from the wind data using finite differences. Then taking 
logarithms of both sides of Eq. (16) and solving for p resulted in the 
following equation. 




)-f A kZp 

A lo 



-f- 





-Z+T-o 



( 20 ) 
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The only variable for a given level is the numerator of the third term 
inside the brackets. 

The final variable to determine for this investigation was the 
friction velocity. The working equation was developed by taking logarithms 
of both sides of Eq . (18) and then taking e to the power of the left side 
and the right side, giving 



Q (21 

This equation looks more difficult to solve than the original equation, 
but is actually more adaptable to programming on our electronic digital 
computer . 

Prior to solving Eqs. (19), (20), and (21) the "constant 11 parameters 

★ 

(von Karman's constant k, adiabatic mixing velocity w , and roughness 

a 

parameter z ) had to be determined. 

★ 

According to [8], w and z are easily obtained by plotting the wind 

a o 

speeds against In z for adiabatic conditions, and drawing a best-fitting 
straight line. The roughness parameter being the intercept on the In z 
axis and the mixing velocity being the slope divided into von Karman's 
constant, or 




t?(U; l- U() 




( 22 ) 



The first problem to be solved was that of determining the time, and 
therefore the wind-sounding, for neutral conditions. All of the layers 
examined up to 110 feet were of depth not exceeding 10 meters. In layers 
of such depth, bounded by thermistors, a dry adiabatic lapse rate 
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corresponds to a vertical temperature difference |aT| ^0.1°C. In such 
shallow layers, a neutral lapse may be regarded as essentially iso- 
thermal. Special emphasis in this regard was given to the temperature dif- 
ference in the layer 0-15 feet above the ground. Runs 16 and 17 corre- 
spond to stabilities nearly neutral for the entire period. 

On clear days, a neutral surface layer is considered [8] to occur 
near sunrise or sunset. This was found to be the case. The actual time 
of isothermalcy usually occurred between half-hourly wind observations near 
sunrise. Hence determining appropriate wind values required the averaging 
of at least two consecutive sets of data centered about the time of iso- 
thermalcy. Actually it was decided to average wind speeds for an hour 
before and after the neutral case in order to reduce round off errors in 
the reported wind values, since these were given only to the nearest mph. 
These averaged wind speeds were then plotted against In z and a straight 
line fitted. The straight line went through or nearly through all the 

plotted points. The values of z q were consistent with those 12 hours 

* 

earlier or later. The adiabatic mixing velocity w still varied from the 

morning to evening adiabatic case. This is, however, to be expected. The 

adiabatic mixing velocity thus obtained was used for six hours preceding 

and six hours after the time of neutral stratification. 

Suggestions have been made [16], [17] that k is a function of stability. 

However, most specialists in the field treat k as a constant independent of 

stability. Values given for it in textbooks and in the literature range 

from 0.38 to 0.45, with most references giving k as 0.40. This is the 

value used in this thesis. The particular constant value used for k will 

2 

not materially affect the results of correlating p with (1+x) . 

t 

It was decided to use the CRA 102-A general purpose automatic digital 
computer [19] rather than a slide rule or hand computer. This decision was 
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based on the large number of computations involved and the additional 
numbers of significant digits possible with the computer. 

There are two possible disadvantages in using such a computer, the 
first being that the data had to be converted from the decimal number 
system to the octal number system. If the data were whole numbers, such 
as the wind observations, this presented little difficulty because programs 
are available for the necessary conversion. However, when involved with 
mixed numbers, as in the temperature observations, it is necessary to con- 
vert the whole number part and the fractional part separately, and then 
add them together. It was possible, however, after a little experience to 
rapidly convert the whole numbers mentally, thus necessitating machine 
conversion of the fractions only. It was found that for a run of 24 hours 
with wind and temperature at four heights and the mean temperatures of the 
three layers (a total of 528 numbers) that the data could fee converted and 
put on IBM cards in less than one hour* 

The second disadvantage was the writing of the programs for the com- 
puter. It is estimated that a total of 20 hours was spent writing the 

2 

necessary programs (Ri, d u , (1+x) , f3 , u^, and a correlation coefficient 

d z 

program), and another 20 hours spent testing these programs on the com- 
puter to make certain they were working properly. 

Once this preliminary work was accomplished the remaining work was 
relatively simple, since the programs were good for any run once a few con- 
stants were changed. An example of the time saved will be illustrated. 

Five friction velocities were computed by hand using logarithms and a slide 
rule, taking nearly one hour. Using the computer, the friction velocities 
for one 24-hour run and three layers (244 values) were computed, typed out 
in both octal and decimal, and IBM cards punched in 15 minutes \ Assuming 
that the time for the hand calculations could be reduced by half, it would 
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still take over 24 hours or about 100 times as long as the machine. 

By having the computer punch IBM cards for all variables it was 
then an easy matter to determine correlation coefficients between the 
various parameters. It took the machine 36 seconds to compute the 
correlation coefficient between two parameters consisting of 48 observa- 
tions each. 

The computer can also plot the information taking a little over a 
minute to plot 48 values of two variables. This was especially useful 
because the relatively small range the various parameters had would make 
it very tedious to plot by hand. 

Once a program has been successfully tested, there is little chance 
of erroneous values. With hand calculations, arithmetical errors may 
occur at any time. 

One last advantage that seems worthy of mentioning is the fact that 
with all the information on cards, it would be easy for another person to 
use in an investigation along similar lines - such as in determining the 
heat flux or the coefficient of eddy viscosity. 
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5. Results 



A discussion of the results of this investigation and a comparison 

of these results with those of other investigators into surface layer 

turbulence will now be made. The numerical values of the parameters 
2 

(1+x) , p , and u* for each run may be found in tables in Appendices I, 

II, III, and IV. It will be seen in these tables that for Runs 16 and 17 

three layers were evaluated - 3 to 12 feet, 12 to 41 feet, and 41 to 

91 feet. Although Deacon's wind profile is usually considered applicable 

only to about 41 feet, a higher layer was analysed in order to determine 

the inconsistencies, if any, from those at lower elevations. Runs 8 and 

10 were used, even though it was only possible to get values for the 

middle layer (due to insufficient data), because of the much greater ranges 

of stability than for Runs 16 and 17. 

2 

The diurnal variation of (1+x) , p , and u^, as well as their 

correlation with one another is illustrated in Fig. 1, where the three 

parameters for Run 8 (12 to 41 feet) are plotted against time. Run 8 was 

taken on 27-28 September with clear skies and light to moderate winds 

2 

prevailing for most of the run. It is seen that (1+x) generally follows 

the diurnal march of temperature, becoming greater than one after sunset 

and less than one after sunrise, whereas f3 has the opposite trend. For 

this particular run f3 was less than one for only a few hours just prior to 

midnight, indicating the layer was unstable for the majority of the run. 

2 

The inverse relationship between p and (1+x) can easily be recognized 
in Figs. 1, 2, and 3. The curves of Figs. 1 and 6 have been smoothed 
somewhat using a smoothing equation 




+ . 25 /? 
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Figure 2. The variation of /? and (l+x) 2 with time of day for 
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Figure 3* Scatter diagram illustrating the variation of ft with (l+x) 2 for 



For a closer examination of the short period fluctuations of j3 , see the 

2 

Appendix. A negative correlation between p and (1+x) is to be expected 
since x is directly related to the Richardson number 



x * 




r; 



The quantities ft and (1+x) in Figs. 2 and 3 are weighted means for the 
entire layer from 3 feet to 91 feet for each observation time of Run 17. 
The equation used to determine this mean being 




and a similar equation for (1+x) . Fig. 3 is a scatter diagram of (3 and 

2 

(1+x) , with the line fitted by eye. The scatter of the points is prob- 
ably mainly the result of the inaccuracy of the original wind measurements 

which were given only to the nearest mph. For Run 17 the linear 

— 2 

correlation coefficient between (3 and (1+x) of -.54 was determined. This 
significant negative correlation is considered to indicate the soundness 
of the Lettau-Mart in theory. 

In 1948, Sutton [18] deduced that the Richardson number must be 
independent of height. However, Deacon, in the same year, found from his 

observations and calculations that Ri varied almost linearly with height. 

2 

Fig. 4 shows (1+x) for each layer plotted against time for Run 16. A 

2 

decrease of (1+x) with height (and thus a decrease of Ri with height) is 

2 

very apparent in this figure. The average value of (1+x) was determined 
for each layer, using the usual equation to find a mean 
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Figure 4. The variation of (1+x) 2 with time, for three layers from Run 16 . 
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The average values were then plotted against height, resulting in Fig. 5. 

2 

The linear decrease of (1+x) with elevation, for this run, substan- 

avg 

tiates Deacon's conclusions* 

In comparing the values of ^3 obtained in this investigation with those 
of previous investigators, they were found to fall within the same range. 
Calder [l] had values of (3 from .75 to 1.4, Deacon (1948) found that (3 
varied from .75 to 1.2, and Priestly and Sheppard [15] determined the 
range of to be from .70 to 1.3. A minimum value of (3 equal to .84 
occurred in this investigation during Run 10, with the maximum value of 
1.4 appearing during Run 8. 



Investigator 


Min {3 


Max (3 


Calder 


.75 


1.4 


Deacon 


.75 


1.2 


Priestly & Sheppard 


.70 


1.3 


Johnston 


.84 


1.4 



Table of ranges of (3 as determined by the 
indicated investigators . 



Table 1 

f . 3 was found to vary considerably from one observation time to the 
next, with a maximum variation of around .2 not uncommon for Runs 8 and 
10. Longley [ll] derived an equation. for (3 that made the usual assumption 
of 2 ^ constant with changes in stability unnecessary. In this equation 
(3 depends only on wind speeds at three levels if the heights are related 
2 

by z * ^2 r = z i r > where r is any ratio. His equation was 
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Using this equation, he computed values of fS from Deacon's data and found 
an average difference from Deacon's fi of from 0.1 to 0.2 depending on the 
heights selected. He also determined (3 from other data for various times 
during both day and night. He found that for daylight hours the range of 
four values of p taken during the same hour might lay between 0.1 and 0.2, 
and for inversion conditions the ranges may even be greater. His final 
conclusions were (1) a possible error of 0.1 or 0.2 does not permit dif- 
ferent types of stability to be clearly distinguished in many cases, and 
(2) it is impossible to calculate ^3 with sufficient accuracy to warrant 
its use in any practical problems. Deacon, in answering Longley's article, 
agreed that ^2 is difficult to measure accurately, and said, 

I have never been so optimistic as to suppose that in practical 
problems it would be profitable to attempt the routine measure- 
ments of (3 in the field. 

It should be remembered that Longley's method of measurement was 

based essentially on curve fitting and had no way of separating out the 

effects of instability. Longley's conclusions may no longer be tenable 

1 

since in the present paper, use of (1+x) accounts for a considerable part 
of the variability of J3 and of u^. 

In perusing the available literature on surface layer turbulence, etc., 
no investigations into the possible variation of {3 with height were found. 

(3 is generally considered constant with height, at least in the lower 
40 feet. However, Deacon (1948) felt there was a possibility that j3 
might vary slightly with elevation, and that if it did, it should slightly 
increase. Fig. 6 indicates a possible variation of /? with height. 
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Figure 6. The variation of ^ with time, for three layers from 
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3 , from Run 16, was averaged for each layer using the same equation 
2 

as that for (1+x) , with the results as below, 

avg 






0.9702 



/? * 1.0040 

/ 2 ^" 



/3 = 1.0271 

A * 6 

It was found that in plotting these values against In z, as in Fig. 7, 
the points could be connected by a straight line. This suggests that (3 
increases with the logarithm of height in the friction layer. The average 
(3 for each layer was also determined for Run 17. Here it was found that 
the average for the upper layer decreased slightly from that for the 
middle layer. However, a straight line connecting /^ V g f° r t ^ ie middle 
and lower layer gave a similar slope as that found for Run 16. This could 
possibly illustrate the fact that Deacon's equation was not applicable to 
as great a height for Run 17 as for Run 16, due to a difference in thick- 
ness of the surface layer, or it could mean that the relationship of 
with height found from Run 16 was only accidental. 

Prior to this paper, there has been no equation for accurately determin- 
ing the friction velocity for non-neutral conditions without measurements of 
the surface drag or short period wind velocity fluctuations. Panofsky [14], 
by making several assumptions, derived an equation for u^ which depends 
only on the wind speed at three elevations. However, upon comparing his 
results for u^ with those obtained by measurements of the short period 
velocity fluctuations, they were found to be less than satisfactory. 

Deacon [3], while making an investigation into the variation of the 
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eddy stress with height, found that values of at 75 feet were much 
more variable than those at five feet. He attributed this mainly to 
accelerations associated with large scale eddies. He also found that 
there seemed to be no systematic variation of u + with height, although 
individual measurements showed large departures. 

Fig. 8 is the friction velocity from Run 16 for the three layers plotted 
against time, as determined by Eq . (18). Here it is seen that the fluctua- 
tions of u^ for the highest layer are generally of greater amplitude than 
for the lower two layers. Although this run seems to indicate a possible 
increase of u* with height, Run 17 indicated a possible decrease with height. 
However, there does not seem to be a systematic increase of u^ between the 
layer 3-12 feet as compared to 12-41 feet. Any possible conclusions that 
can be made here regarding the variation of u^ in the layer 41-91 feet are 
tentative, inasmuch as the Deacon profile may not be applicable to these 
elevations. However there is no clear-cut increase of u * between the layer 
12-41 feet and 41-91 feet. 
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Figure 8. The variation of the friction velocity with time, for three layers from 



6. Summary and Conclusions 

Martin’s combination of Lettau's vertical nixing velocity with 
Deacon's wind profile equation gives values of p and which compare 

favorably with those of other investigators. However, the equations con- 

2 

taining the term (1+x)^ will not apply for x< -1, which may have the 
significance of placing an upper limit on the depth of the surface layer. 
During this investigation x was never found to be less than minus one. 

The equation for the friction velocity, Eq . (18), may become an 
accurate method for determining u^ if it becomes possible to evaluate (3 
with sufficient accuracy. Additional work is needed to compare u*, as 
determined from Eq . (18), with that obtained from surface drag measurements 
or short period wind velocity fluctuations, before an accurate evaluation 
of the equation can be made. 

By comparing a large number of cases it was found possible to use data 
that was otherwise undesirable because of the inaccuracies of the wind 
measurements. However, to determine the various parameters for any one 
observation, it is felt that wind measurements of much greater accuracy are 
needed . 

One of the most important results of this investigation is the determin- 
ing of a possible variation of p with height and the correlation of ft with 
2 

(1+x) . However, further work with a larger number of layers is needed to 
verify this result. 
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Appendix I 
Run 8 



time 


(l+xf 


/? 




1300 


.528 


/ 

1.133 


.899 


1330 


.496 


1.226 


.583 


1400 


.496 


1.226 


.583 


1430 


.497 


1 . 141 


.919 


1500 


.598 


1.168 


.661 


1530 


.633 


1.161 


.648 


1600 


.708 


k. 147 


.623 


1630 


.708 


1.119 


.721 


1700 


.827 


1.163 


.489 


1730 


.911 


1.201 


.362 


1800 


1.000 


1.275 


.223 


1830 


1.399 


1.098 


.408 


1900 


1.684 


1.076 


.383 


1930 


2.196 


1.007 


.422 


2000 


2.336 


.972 


.478 


2030 


' 2.269 


.976 


.483 


2100 


2.487 


.964 


.468 


2130 


2.208 


1.007 


.421 


2200 


1.879 


.976 


.481 


2230 


1.945 


1.022 


.440 


2300 


2.078 


1.014 


.430 


2330 


1.758 


1.007 


.527 


0000 


1.408 


1.062 


.492 


0030 


1.638 


1.043 


.467 


0100 


1.464 


1.030 


.562 


0130 


1 .464 


1.093 


.402 
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Run 8 ( cont .) 



time 


0+xf 


/9 


u* 


0200 


1.408 


1.098 


.407 


0230 


1.522 


1.088 


.396 


0300 


1.409 


1.097 


.407 


0330 


1.145 


1.173 


.335 


0400 


1.196 


1.118 


.431 


0430 


1.196 


1.118 


.431 


0500 


1.249 


1.077 


.512 


0530 


1.302 


1.157 


.320 


0600 


1.146 


1.123 


.437 


0630 


1.125 


1.067 


.373 


0700 


.882 


1.147 


.311 


0730 


.670 


1.275 


.257 


0800 


.409 


1.327 


.257 


0830 


.492 


1.303 


.178 


0900 


.216 


1.405 


.320 


0930 


.277 


1.375 


.294 


1000 


.497 


1.217 


.379 


1030 


.221 


1.403 


.318 


1100 


.347 


1.212 


.560 


1130 


.223 


1.266 


.652 


1200 


.224 


1.266 


.652 


1230 


.089 


1.379 


.893 
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Appendix II 



Run 10 



time 


mr 


-A- 




1630 


.459 


1.103 


.627 


1700 


.575 


1.081 


.565 


1730 


.574 


1.081 


.566 


1800 


2.716 


.873 


.368 


1830 


1.112 


.956 


.575 


1900 


2.217 


.854 


.499 


1930 


2.470 


.844 


.472 


2000 


1.552 


.989 


.344 


2030 


.032 


1.150 


1.610 


2100 


.180 


1.124 


1.430 


2130 


2.233 


.891 


.406 


2200 


2.488 


.881 


.384 


2230 


1.557 


.989 


.343 


2300 


2.493 


.945 


.271 


2330 


1.359 


.937 


.520 


0000 


1.173 


.951 


.560 


0030 


1.000 


.966 


.606 


0100 


1.361 


.937 


.520 


0130 


1.174 


1.015 


.396 


0200 


1.440 


.918 


.704 


0230 


1.563 


.910 


.675 


0300 


1.441 


.918 


.704 


0330 


1.324 


.963 


.599 


0400 


1.566 


.948 


.551 
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Run 10 (contd.) 



time 


il±2L t 


0430 


1.324 


0500 


1.000 


0530 


.559 


0600 


.486 


0630 


1 .444 


0700 


1.325 


0730 


.111 


0800 


.423 


0830 


.423 


0900 


.492 


0930 


.642 


1000 


.165 


1030 


.257 


1100 


.368 


1130 


.498 


1200 


.048 



/? 




/ 

.963 


.599 


.952 


.845 


1.043 


.922 


1.056 


.989 


.918 


.703 


.926 


.734 


1.019 


.812 


1.133 


.750 


1.069 


1.060 


1.119 


.695 


1.095 


.608 


1.220 


1.200 


1.201 


.856 


1.082 


1.137 


1.054 


.977 


1.089 


4.596 



* 
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APPENDIX III 







X 


o 


CM 


ON 


oo 


r- 


ON 


CM 


CM 


CM 


O 


m 






O 


— 1 


o 


co 


O' 


X 


oo 




r—4 


oo 


CO 


00 




1 ^ 


O 


o 


crx 


O' 


O' 


O' 


ON 


o 


o 


• o 


o 


O' 






. — i 


•— < 














pH 




r—4 






CM 




























✓— s 


CM 


CO 




O' 


CO 


r—4 


o 


<f 




00 


o 


X 




X 


00 


CM 


00 


m 


m 


CM 


—4 


X 


<r 


i 


CO 


r- 




*t 




00 


00 


oo 


O' 


On 


o 


X 


X 


r>. 


oo 


r-^ 




*— 4 




















• 


• 






























































X 




























Gl 


co 


oo 


CM 




<3 


CM 


m 


CO 


ao 


r— 4 


CO 


r—4 




s 


CO 


f-H 


O' 


X 


CM 


CO 


00 


X 


<3 


m 


oo 


m 






X 


m 


<3 


X 


CM 




o 


*>0 


<3 


■ oo 


CM 


00 




* 


• 


















• 


• 


• 




3 


f-H 




rH 




f—4 


-4 


—i 


— 1 


pH 


p— 4 


r—4 


r—4 


—1 




























O' 




























1 




CO 


X 


•<r 


m 


<3 


Mf 


—4 


X 


o 


oo 


CO 


oo 






O 


O 


O' 


O 


i—4 


X 


X 


Mf 


CO 


oo 


00 


oo 


^Hl 


\ 


^ o 


o 


O' 


O 


o 


O 


On 


O 


O 


l o 


o 


o 
























• 




• 






r—H 






—4 


. — 1 












i—4 




u 




























<u 




























>> 






























CM 


























X 




CM 


CM 




m 




X 


O 


X 


—4 


CO 


m 


CO 




X 


r-v 


CO 


<r 


m 




CM 


X 


CM 


oo 


• X 


o 


X 




+ 


X 




r^* 


X 


oo 


00 


00 


X 




X 




X 




. — i 




















. 


* 


. 


1 


























































X 


Mi- 


m 


00 


00 


m 


O' 


X 


CM 


m 


o 


oo 


CM 




a 


en 


CO 


X 


o 


CO 


CO 


—4 


CO 


X 


X 


CM 


X 




E 


r-4 


O 


o 


CO 


CM 


CM 


On 


r—4 


CO 


i m 


Ml- 


m 




'w' 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 


• 




* 


r— 4 




i—4 


«— 4 


r— 1 


r—4 




"4 


r-4 


r—4 


r—4 


—4 




3 


























1— < 
























































i 




o 


X 


in 


oo 


O' 


O' 


CM 


X 


o 


r—4 


X 


—4 


- 


<y> 


o 




On 


CO 


X 


X 


mT 


o 


o 




o 


P"- 


CM 

y 4 


X 


. °. 


o 


On 


O' 


O' 


O 


o 


o 


o 


1 O 


o 


o 






r— ^ 


r—4 










f-H 


r— 4 






































<3 
























































CS 


CM 


























X 


/-v 




H 


O 


r^ 




O' 


o 


<r 


CO 


o 


ON 


00 




X 


o 


CM 


X 




m 


<3 


X 




r—4 


o 


m 


o 




+ 


O 




o 


i—4 


o 


o 


O' 


o 


o 


1 00 


oo 


OO 
























• 


• 


• 










1 1 1 


p—4 


H 


















y-*S 




























X 




























CL 


<* 




oo 


m 


CM 


CM 


CM 


^-4 


m 


o 


1—i 


p^ 




£ 


uO 


O 


oo 




CO 


ON 


o 




m 


m 


ro 


00 




w 


o 


CM 


r^ 


oo 


O 


r—4 


r-H 


o 


CM 


1 CM 


m 


o 





3 






r—4 




p—4 


-4 


p—4 


p—4 


^-4 




f-H 


pH 


f— H 


CM 






























i —4 




CO 


X 


CM 


m 


00 


CO 


CO 


o 


m 




3 


X 


O 


1 




pH 


p^ 


P^ 


X 


o 


r- 


p^ 


r—4 


oo 




CO 


CM 


CO 


O 


o 


o 


o 


o 


o 


o 


O 


o 


l 


o 


ON 


o 




—4 






r—4 


r-4 






—4 










p—4 


L4 






























o 




























































3 


CM 




























X 




CM 


X 


CM 


CO 


mT 








m 




?~4 


i—4 


vt 




X 


O 


o 


MJ- 


X 


CO 


ro 


CO 


r—4 


r—4 




CM 


CM 


1—1 




+ 


O 


o 


o 


o 


O 


O 


o 


o 


O 


l 


o 


O 


O 






• 




















• 


• 


• 








1—4 






r^ 


p— H 


r—4 


r-H 


—4 




r-4 


r—4 


pH 




a> 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


o 


O 




E 


CO 


o 


CO 


o 


CO 


o 


CO 


o 


CO 


o 


CO 


o 


CO 




•r>4 


X 




p^ 


OO 


oo 


o 


o 


o 


O 


—4 




CM 


CM 




4-4 


r—4 


p—4 


p—4 


p—4 


— < 


—4 


1—4 


CM 


CM 


CM 


CM 


CM 


CM 






38 



h 




<* 


00 








CO 


00 


to 


NO 


<* 




rH 


^4 


rH 




O 




co 


o> 


CM 


r— * 


00 


o 


. — I 


CM 


i— g 


oo 


oo 


Mf 






O 


on 


O 


CO 


o 


o 


CO 


o 


o 


O 


o 


o 


O 


O 










r— ^ 




— 


- 








- 


rH 


- 


rH 


r-H 


leg 
































/-N 




ro 


— * 


00 


<* 




X 




CO 




1/1 


uo 


oo 


f — i 




X 




CM 


00 






r— t 


0 o 


oo 


CM 


o 


. 1 


r^ 


00 


vO 


ON 


t 




r^* 






00 


00 


00 


O'- 


00 


00 


CO 






r*> 








* 




























































✓“N 
































X 




o 


v£> 




Mf 


ON 


00 


^H 


CT\ 


00 


UO 


oo 


oo 


ON 


X 


a 




r— 1 


<* 


on 


u0 


00 


Mf 


uo 




oo 


. — 4 


oo 


oo 


rH 


O 


E 




o 


ON 


vO 




u0 


<f 


00 


uO 


uo 


uo 


vO 


vO 




rv 


V ' 




• 




























4 




CM 


rH 




rH 




rH 


^-4 


r— 4 




rH 




_ 


rH 


—4 


3 







































uo 


— i 


vf 


CM 


uo 


O 


oo 


o 


uo 


<fr 


rH 


rH 


nO 


UO 






ON 




Mf 


ON 


oo 


oo 


00 


o 


oo 


oo 






Mf 


Mf 


P 




on 


ON 


O 


ON 


o 


o 


ON 


o 


o 


O 


O 


o 


O 


o 


0) 






• 


* 


























so 

c 










f a 4 






*"* 




-* 




i-H 




r— < 


—4 






CM 


































/-V 




CM 


. — I 


-H 








CO 


ON 


00 




so 


X 


00 


o 




X 




O 


uo 


UO 


oo 


CM 


00 


vO 


uO 


CM 






r^ 


oo 






+ 




vO 


vO 


vO 


00 






vO 








vO 


v£> 


nO 


X 









































y—s 


































JC 


<* 


o 


OO 


CM 




vO 


CM 




CM 


CM 


UO 


oo 


CO 


CM 






a 


uo 




00 


o 


o 


CM 


r^ 


^4 


rH 


CM 


00 


o 


fH 


CO 






E 


CM 


oo 


On 


vf 


CM 


CM 


oo 


CM 




CM 


uo 


CM 


CM 


o 


✓— 0. 




'w' 




























• 


• 




4 




^4 




^4 


^—4 


, — 4 


, — 4 


r-4 


, — 4 


, — 4 


f— 1 


r— t 


—A 


r— 1 


J-4 




3 






























c 


- 
































o 


■ ^ 
































o 


vf 
































v_^ 


» 




uO 


oo 


uO 


CM 


r* 


*— 4 


00 


ON 




o 


vf 


00 


ON 


Mf 




- 


to 


CM 


00 


uo 


ON 


i — i 


CM 


00 


—1 


ON 


CM 




i—4 


r-4 


vO 


vO 


CM 


X 


O 


On 


O 


ON 


O 


O 


ON 


o 


ON 


O 


ON 


o 


O 


O 


f- 4 


r— 4 






































, — 4 




i—4 




~—4 


, — 1 




, — 4 




, — j 




r—i 


^_4 


r—< 


c 


u 
































3 


0) 
































c£ 


So 


































CC 


CM 
































*-3 


x-s 


UO 


OO 


O'* 


rH 




oo 


i—4 


o 


00 


o 


CM 


ON 


CM 


<t 






X 


OO 


OO 


o 


ON 


o 




oo 


ON 


o- 


00 




On 


ON 


On 






+ 


00 


ON 


ON 


00 


ON 


oo 


ON 


00 


00 


00 


oo 


00 


00 


00 









































































X 


CM 


uo 


—A 


vD 


vO 


00 


uo 


uo 


o 


uo 


r- 


o 


© 


o 




a 


vt 


uo 


CO 




vf 


ON 


uo 


UO 


uo 


uo 


00 


ON 


uo 


ON 




Q 


CM 


CM 


uo 


CM 


CM 


co 


CM 


CM 


CM 


CM 


o 


o 


CM 


uo 




S-/ 
































4= 




r—4 


f—i 


rH 


1—4 


r-»4 


r— 4 


*—A 


rH 


rH 


rH 


rH 


rH 


rH 




J 


































co 


UO 


\D 


3 




00 


UO 


uo 


t 


uo 


o 


rH 


Mf 


uo 


CM 




00 


00 


CO 


uO 


00 


00 


oo 


CM 


CM 


00 




l 




ON 


ON 


ON 


ON 


ON 


ON 


On 


ON 


ON 


ON 


o 

i—4 


O 


ON 


ON 


co 
































u 
































0) 
































So 


CM 


vD 


<t 


r-i 


ON 


ON 


r-H 


<t 


Mf 


rH 


<r 






rH 


O 


<TJ 




co 


r^4 


CM 


CM 


CM 


CM 


H 


rH 


CM 


rH 


rH 


o 


CM 


co 




X 


O 


o 


O 


o 


o 


O 


o 


o 


O 


o 


o 


o 


O 


o 




-+ 


• 






























—4 


H 


—4 


i—4 


1—4 


—4 


r— 4 


i—4 


— g 


i—4 


rH 


r*H 


rH 


i—4 


—4 




































V 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 




B 


o 


CO 


o 


CO 


o 


CO 


o 


co 


o 


co 


o 


co 


o 


ro 






CO 


CO 


o 


o 


i—4 


—A 


CM 


CM 


ro 


ro 






uo 


uo 




U 


CM 


CM 


o 


o 


o 


o 


O 


o 


o 


O 


o 


o 


o 


O 



39 




40 






810 ’ 



APPENDIX IV 







r>. 


m 


Mf 




CNI 


r~i 


o 


vO 


oo 


Os 


CNI 


o 


CO 






os 


sO 


sD 


vO 


OS 


o 


Mf 




sO 


o 


m 


F— *4 


<r 




ON 


ON 


ON 


o> 


Os 


o 


Os 


os 


Os 


o 


crs 


o 


Os 
















-• 








i— 4 




r— l 




CM 


































co 


Mf 


00 


CM 


CO 


o 


00 


•— 1 


sO 




O 


—4 


sO 


X 




os 


CNI 


s£> 




sO 


vO 


00 




O' 


Mf 


' Mf 


OS 


in 


+ 




vO 


00 


00 


O 


Os 


Os 


Os 


Os 


00 


Os 


Os 


Os 


Os 


. — t 








































— * 





















x: 





cl 


CNI 


CNI 




00 


r>. 


CM 


—4 


m 


r-. 


CM 




m 


r-s. 




to 


ON 


i— 1 


CO 


CM 


o 


Mf 


CM 


CO 


so 


Mf 


r>. 


oo 


s£) 






CM 


CNI 


—4 


Os 


o 


Os 


O 


p'4 


•—1 


Os 


o 


00 


O 


- 


* 




























»— i 




«— I 


•—4 


r— ^ 




r-i 




» —4 


r— * 


p^ 








f— 4 


Os 

I 


































Mf 


Os 


in 


i—4 


sO 


co 


M0 


m 


v£> 


CO 


o 


00 


co 


Mf 






CM 


c- 


co 


Os 




co 


r- 


m 


>3- 


m 


CO 


CO 




X 


o 


Os 


Os 


Os 


Os 


o 


os 


Os 


Os 


o 


os 


o 


OS 


U 




N 


























<U 




< 










i—i 








—4 




1—4 




Pn 






























tc 
































CNI 


































CM 


i— 1 


o 


Os 


00 


oo 


o 


O 


oo 


oo 


CO 


o 




X 


Os 




Mf 


o 


*— 4 


r— 4 


os 


Mi- 


>$ 


*— < 


f— 4 


00 


in 




+ 


sO 


00 


00 


o 


Os 


Os 


OS 


oo 


00 


Os 


Os 


OS 


O' 




.—4 




























































/-\ 






























x: 






























O. 


sO 


CM 


sO 


m 


00 


00 


m 


m 


00 


CM 


Mf 


CO 


Mf 




6 


m 


CM 




r— 4 


o 


CM 


c- 


CM 


x> 


m 




CM 


00 




N-/ 


m 


CM 


O 


Os 


Os 


Os 


Os 


Os 


Os 


os 


Os 


Os 


Os 




■K 


































r— 1 


r-^ 






















X 






























Mf 

i 




































o 


sO 


sO 


m 


r>. 


r-4 




o 




oo 


00 


Os 


CNJ 


co 


oo 


Mf 


m 


Mf 


oo 




m 


oo 


Os 


sO 


sO 


•—4 


sO 


. — 1 


X 


os 


o 


os 


Os 


os 


Os 


Os 


os 


Os 


Os 


Os 


o 


Os 


iH 






r— 4 




















H 




<u 






























Po 






























dJ 






























J 


CM 






























/N 


Mf 


co 




CO 


OS 


Os 


Os 


o 






Mf 


—4 


m 




X 


Os 


Mf 


sO 


CM 


?— 4 


o 


m 


o 


m 




m 


Os 


Mf 




-f 


m 


c- 


oo 


O 


o 


o 


os 


o 


Os 


Os 


os 


Os 


Os 




— 1 






























w 








f— ^ 


r— 4 


p—4 


















/N 






























XL 






























CL 


m 


in 




co 


vf 


00 


m 


Mf 


Mf 


i—i 


CM 


CM 


oo 




E 


sO 


x> 




, — i 


Mi- 


«— 4 


Os 


CM 


Os 


O 


in 


Mf 


Mf 




V-/ 


CO 


CO 


co 


CM 


os 


CM 


o 


CM 


o 


>—4 


Mf 




Mf 




* 






























3 






r-4 








p^ 




1—1 


—* 


f— 4 


p—4 




CM 




























































i 




00 


oo 


m 


o 


Os 


•—4 


CO 


CM 


CO 


Mf 


r-^ 


Os 


o 


- 




CNI 


CM 


CM 


Mf 


00 


Mi- 


sO 


Mf 


vO 


sO 


«— 4 


o 


p^ 


co 


Os 


Os 


os 


Os 


Os 


os 


Os 


Os 


Os 


Os 


Os 


Os 


Os 


u 






























(V 






























rO 






























J 


CM 
































00 


00 


'O 


CM 


p*4 


Mi- 


CO 


Mf 


Mf 


CM 


00 


CO 


Mf 






00 


00 


CM 


Mf 


co 


en 


CM 


CM 


CM 


r-^ 


— 4 


CO 


CM 




+ 


OS 


Os 


O 


o 


o 


o 


O 


O 


O 


o 


o 


o 


O 




•—4 








































p'4 




rW 


r— 1 














<u 


o 


O 


o 


o 


o 


o 


O 


o 


O 


o 


o 


o 


o 




e 


CO 


o 


CO 


o 


cn 


o 


CO 


o 


CO 


o 


CO 


o 


CO 






v£> 


Cv 




oo 


00 


OS 


Os 


o 


o 


1—1 


1—4 


CM 


CM 




XJ 


i — i 




— 1 




— > 






CM 


CM 


CM 


CM 


CM 


CM 



4 I 



2300 1.028 .974 1.162 .925 1.002 1.125 1.013 .998 1.031 -986 .997 



17 (Cont. 







ON 


CNI 


00 


r*— 4 


00 


M0 


o 


CN 


r>- 


o 


m 


m 


m 


NO 


Mf 






U0 


00 


MO 


oo 


00 


ON 


r» 


r>- 




00 


00 




CM 


co 


co 




CN 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


o 


O 


o 




































CN 






































ON 


m 


— * 


co 


ON 




r>. 


co 


m 




O 


MO 


r>- 


CN 


—4 


X 




r>- 


O 


CM 


■ 1 -4 


CN 


ON 


r>> 




m3* 


M* 


Mf 


r^- 


00 


*— 4 


00 


■f 




o 


*— i 




*— *4 


o 


ON 


O 






o 


o 


o 


ON 


o 


ON 


*— 4 
































. 


w 




— 1 


»— 4 


— t 


*— 4 


*— 4 




*—■4 


•— 4 




—1 


r—4 






r— 4 







^-s 








































































a 




CM 


r^ 


co 




ON 


r^- 


<r 




CN 


i— t 


MO 


o 


M0 


Mf 


in 




E 




O 


r— 4 


Mf 


Ml* 


r—4 


O 


CN 


r—4 


00 


co 


00 


ON 


in 


NO 


00 




' — ✓ 




O 


ON 


ON 


ON 


o 


r—4 


o 


ON 


oo 


o 


o 


ON 


o 


ON 


ON 


- 


* 




• 






























r—4 






f— 4 








r—4 




f— 4 






f—t 


r—4 




r—4 






ON 

1 




































r— 4 






in 


ON 


o 


r—4 




Mf 


00 


o 


r^- 


00 


CN 


m 


o 


in 


in 


Mf 






r—4 


m 


Mf 


MO 


ON 


O 


NO 


M0 


oo 


ON 


r* 


M0 


o 


Mf 


On 


u 


ON 


ON 


ON 


ON 


ON 


o 


ON 


ON 


ON 


ON 


ON 


ON 


o 


o 


ON 


<U 




































>N 




































CO 




































-4 


CN 




































/N 




NO 


m 


m 


CM 


NO 


MO 


NO 


Mf 


Mf 


co 


ON 


00 




co 


nO 




X 




r>- 


o 


o 


NO 


CN 


co 


NO 


O 


o 


r—4 


ON 


o 


00 


r—4 


MO 




+ 




•—4 


CN 


CM 


r—4 


O 


ON 


o 


CM 


CN 


o 


ON 


r—4 


ON 


o 


O 








*“ i 


r-t 


«— t 


r—4 


r-H 




l-H 




r—4 


r—4 








*— 4 






u 

<u 

p> 

c5 



/^V 

o- 


Mf 


Mf 


in 


O 






oo 


Mf 


Mf 


00 


r—4 


m3- 


o 


O 


£ 


r—4 




CM 


00 


O 


in 


Ml* 


NO 


r^ 




co 




r>- 


00 


s— ' 




O 


o 


ON 


o 


ON 


ON 


ON 


ON 


ON 


ON 


o 


ON 


ON 


0 




r—4 


r—4 




r—4 














r—4 











cn 


o 


nO 


co 


co 


ON 


nO 


CN 


00 


MO 


ON 


CM 


CO 


ON 






co 


cn 


CN 


CM 


ON 


00 


MO 


CN 


NO 


nO 


r-4 


r^- 


MO 


MO 


o 


o 


o 


O 


ON 


ON 


ON 


O 


ON 


On 


o 


ON 


ON 


ON 






r—4 


—4 


r—4 


r—4 








r—4 






f-H 









r>- 

oo 

o 



ON 

n 

ON 




CM 



co 

U 
V 
p N 
03 
-1 





CM 


in 


n 


CN 


r* 


in 


MO 


ON 


r—4 


On 


co 


CN 


co 


n- 


Cn. 


M0 


o 


CM 


O 


co 


MO 


O 


ON 


o 


ON 


ON 


CN 


n 


n 


Mf 




co 


M3- 


m 


in 


r—4 


o 


CN 


co 


CN 


CN 


00 


r—4 




. 


^4 


H 




r—4 


r—4 


*— 4 


*— 4 


*— 4 


r-^ 




r—4 




^-4 


— * 


— * 





Mf 


MO 


Mf 


r—4 


n 


r- 4 


co 


co 


co 


CM 


ON 


co 


ro 


CM 


n- 


n 


co 


in 


CM 


r* 


O 


in 


in 


n 


n 


co 


r>- 


n* 




ON 


ON 


ON 


ON 


ON 


ON 


o 


ON 


ON 


ON 


ON 


o 


ON 


ON 


On 



CM 

/"-V 


ON 


o 


MO 


Mf 


00 


co 


00 


00 


M0 


00 


00 


00 


Mf 


00 


in 


X 


CM 


CM 


CM 


CM 


CM 


CN 


CN 


CM 


CM 


CM 


CO 


CM 


Mf 


CO 


uo 


+ 


O 


O 


O 


O 


O 


O 


O 


o 


O 


o 


O 


O 


o 


o 


o 


»— 4 






























• 






r—4 


r—4 




r—4 




r—4 


r—4 


t— *4 


r—4 


rr4- 




*— 4 


r—4 


r—4 



o 


O 


o 


O 


o 


o 


o 


o 


8 


o 


o 


O 


O 


o 


O 


o 


E 


CO 


o 


CO 


o 


CO 


o 


CO 


CO 


o 


CO 


O 


CO 


o 


CO 


•t-4 


CO 


O 


o 


1—4 


r—4 


CM 


CM 


CO 


CO 


m3- 


M3- 


in 


n 


MO 


MO 


u 


CM 


O 


o 


o 


o 


o 


o 


o 


O 


O 


O 


o 


o 


O 


O 



42 



1 




17 (cont. 






<fr 


on 


<5- 




00 


ON 




rH 


ON 


CM 


<r 






ON 


r-H 


CM 


un 


p* 


00 


ON 


ON 


00 




H 


co 


o 


o 


co 


r-H 


O 


o 


ON 


ON 


00 




ON 


ON 


o 


o 


o 


r— 4 


o 


r-H 


-H 


— A 
















r— 4 




r-H 


r-H 



rH 

CM 

o 



cm 


X 


CM 




o 


co 


CM 


Mf 


ON 


O 


m 


ON 


ON 


un 


o 


r-H 


X 


O 


ON 


ON 


co 


00 


o 




m 


r-H 




X 


m 








-t 


o 


ON 


ON 


o 


ON 


o 


ON 


ON 


o 


r-H 


00 


n* 


oo 


00 


p^ 


» ■■ J 


































fH 






r-H 




— < 






rH 


r-H 













X 





a. 






ON 


CM 


CNJ 


ON 


co 


CO 




p* 


m 


ON 


o 


CO 


o 




E 


CM 


Mf 


o 


m 


o 


r-H 


x 


o 


ON 


oo 


m 


Mt 


X 


o 


r-H 




NH 


00 


r-H 


o 


o 


r-H 


CM 


m 


CM 


ON 




r-H 


m 




CO 


r- 


— 


* 


• 






























r-H 


3 




r-H 


r-H 


r-H 


r-H 


r-H 


. — i 


r-H 






r-H 


r-H 


r-H 


r-H 


r-H 



on 

i 



r-H 




CO 


x 


r*- 


r-H 


CO 


-H 


00 


oo 




o 


oo 


00 


Mt 


m 


Mt 




(V, 


CO 


n- 


ON 


m 


p** 


00 




m 


>d- 


CO 


m 


CM 


CM 


r-H 


o 




X 


o 


ON 


ON 


o 


ON 


ON 


ON 


ON 


O 


o 


o 


O 


O 


o 


r-H 


u 


> 


• 


• 




























a> 




r-H 






r-H 










^H 


r-H 


rH 


r-H 


r-H 


f-H 


r-H 




































CO 


































X 


CM 




































m 


r-H 


X) 


ON 




m 


o 


ON 


r^ 




r-H 


o- 


f-H 


o 


m 




X 


o 


r-H 


r-H 


00 


m 


<5- 


Mf 


00 


m 


NO 


ON 


ON 


MT 




CM 




+ 


CM 


o 


r-H 


ON 


o 


crs 


r- 


ON 


o 


CO 


00 


NO 


r- 


00 


NO 




f-H 


# 


o 


• 
























• 




'w' 




r-H 


r-H 




rH 








f— H 


r-H 













X 

a. 


n- 


o 


r-H 


nO 


CM 


00 


ON 


ON 


CO 




rH 


Mf 


m 


ON 


CO 


B 


p^ 


m 


00 


o 


r-H 


oo 


Mf 


nO 


CO 


nO 


ON 


ON 


CO 


CO 


CO 




o 


o 


r-H 


CM 


CO 


CO 


CO 


f-H 


CO 


CO 


CO 


CO 


CM 


CO 


M* 


* 
































D 


-h 


r-H 


r-H 


r-H 


r-H 


r-H 


r-H 


f-H 


rH 


r-H 


rH 


r-H 


r-H 


f-H 


r-H 



c 

3 

OZ 



I 

N 



U 

<D 

>y 

05 

XJ 




iO 


00 


CO 


r^ 




o 


X 


m 


f-H 


*n 


ON 


CJN 


r-H 


^H 


oo 


r- 


CM 


00 


o 


rH 


m 


r^- 


nO 


nO 


ON 


rH 


p*. 




o 


ON 


ON 


O 


ON 


o 


o 


o 


ON 


ON 


ON 


ON 


o 


ON 


o 


rH 


ON 




rH 




r-H 


^H 


r-H 










r-J 




r-H 


r-H 





CM 


m 






X 


o 




NO 


ON 


f-H 


CM 


o 


rH 


m 


ON 


CM 




n* 


m 


m 


<N 


m 


ON 


co 


NO 


m 


CM 


o 


o 


ON 


o 


00 


+ 


ON 


o 


ON 


ON 


oo 


p-* 


oo 


ON 


00 


00 


00 


oo 


CO 


00 




r-H 




































rH 


/ 



























x 





a 


CO 


CO 


o 


o 


r-H 




CM 


ON 




00 


CM 


n- 


o 


00 


00 




B 


00 




o 


H 


nO 




r-H 


o 


ON 


00 


CM 


ON 


X 


o 


rH 




V-/ 

* 

0 


CM 


CO 


CM 


CM 


<r 


co 


m 


un 


CO 


CO 


un 


CO 




<r 


un 


- 




rH 


f-H 




^H 


r-H 


rH 


rH 


r-H 


rH 


pH 


rH 


r-H 


rH 


^H 


CM 


































rH 


































- 




rH 


ON 


o 


CM 


r-H 


<t 


r^ 


P* 


o 


X 


ON 


r^ 


CM 


00 


X 


CO 


Co 


un 


un 


oo 


oo 


•<r 


X 




Mf 


NO 


X 


Mf 


X 


CM 


X 






X 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


U 


N 


^ • 






























<D 


































Pn 


































Q 


































X 


CM 


































/-N 




M* 


r>- 


O 


ON 


MT 


• o 


Mf 


o 


00 


p^ 




Mf 


•H 


CM 




X 


un 


CO 


CO 


CM 


Mf 


00 


00 


00 


X) 


X 


X 


m 


X 


mT 






+ 


o 


o 


o 


O 


O 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 




r-H 


• 


o 
































rH 


f-H 


r-H 


^■H 


r-H 























0) 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


£ 


o 


CO 


o 


CO 


o 


CO 


o 


CO 


o 


CO 


o 


CO 


o 


CO 


o 


•rH 




r%. 


oo 


oo 


ON 


ON 


o 


o 


r-H 


f-H 


CM 


CM 


CO 


CO 




Hi 


o 


O 


o 


o 


o 


o 


^H 


rH 


—< 


r-H 


r-H 


r-H 


r-H 


rH 


rH 



43 



1430 .930 .969 1.417 .880 1.011 1.267 .522 1.094 1.873 .682 1.054 




CM 


nO 


on 


— 4 


vO 


<r 


1-0 


m 




MD 


in 


CM 


CM 




00 


ON 


co 


f — 1 


— 1 


CM 


CNJ 




00 


vD 




m 


in 


<r 


ON 


ro 




co 


MD 




C 


o 


O 


o 


o 


on 


on 


ON 


o 


ON 


ON 


ON 


o 


o 


on 


ON 


ON 




r>. 


•— i 


in 




o 


vO 


CM 


o 






vO 


MD 


CM 


CO 


MD 




ON 


ON 




o* 


O 


"O 


ro 


00 


vO 


oo 


00 


co 




in 


<——4 


• — 1 




00 




x> 




ON 


ON 


O 




-O 


co 


in 


— 4 


-< 


CO 


CM 


CO 


<r 


m 














— • 


i— 4 




_< 




— f 






-4 


—4 


•—4 



c 

© 

u 



c 

p 

oc 



u 

Q) 

C 3 

-3 






l-i 

<L> 

Pc 



CM 



co 

U 

© 

?o 

<c 







































-C 




































a 


r^v 


n. 




MT 




CO 


un 


MO 


1— ( 


ON 


ON 


00 


ON 


MD 


o 


00 


r-H 


g 


i—i 


CM 


Mf 


cO 


CO 


r*A 


nO 


vO 


o 


o 




ON 


CM 


MD 


CM 


CM 


o- 




m 


CM 


m 


CM 


CM 


00 


vO 


un 


un 




MD 


un 


un 


un 


un 


o 


co 


p 


— 1 


- 




4 






























r>* 


•— < 


00 


CM 


CM 


*—4 


VO 




un 


o 


MD 


ON 




<3“ 




CO 


Mf 


CVx 


ON 


ON 


CM 


— • 


r-*4 


00 


vO 


CM 


Mf 


CO 


ON 


o 


00 


CM 


»—4 


CO 




\ 


ON 


ON 


o 


O 


o 


ON 


ON 


ON 


ON 


ON 


00 


o 


00 


ON 


ON 


On 


ON 




N* ' 




-< 


i—4 


-« 














- 












CM 






































CM 


CM 


in 


CM 


CM 




—4 


CO 


o 


ON 


co 


<r 


<r 


un 


CM 


ON 


ON 


X 


<r 


ON 


ON 


ON 


ON 




CM 


CM 




00 


oo 


r*. 


CO 


CM 


4"*. 


00 


un 


■t 


r^. 


00 


MO 


00 


00 


ON 


CM 


un 


MD 


o 


CO 


CM 


r^. 


un 


NO 


ON 


«— 4 


r— < 






























• 




. 


'w' 


















— < 


CM 




*— 4 






r^4 


— 1 


CM 


/*\ 




































X 

a 


ON 


r-"4 


un 


in 


un 


CM 


ON 


ON 


un 


CO 


co 


O 


un 


un 


CM 


Mf 

ON 


CO 

co 


g 


o 


m 


ON 


r-4 


O 


00 


00 


00 


O 


un 


ON 


CM 


i—i 


O 


CO 




(—1 


w 


CO 


m 


CM 


CM 


i—4 


MD 


MD 


MO 


00 


ON 


o 


ON 


o 




o 




NJ 

o 


* 

p 


»— 1 


- 1 


“< 


• ■ 4 


r— 4 












—4 






i—4 


—4 


—4 


—4 




<T 




un 


o 


CM 


un 


Mi- 




MD 


ON 


i—4 


MD 


MD 


i—4 


co 


r — 4 

co 


CO 


Oo 




CM 


>* 


N* 


CO 


ON 


vO 


vO 


r^ 


00 


co 


00 


CM 


o 




o 




X 


^ o 


O 


o 


O 


O 


ON 


ON 


ON 


ON 


ON 


o 


ON 


O 


o 


ON 




ON 
































* 


r-L 


* 










»— 4 


f— 4 












•—4 




i—< 


r^4 








CM 


CM 


MD 


CM 


ON 


CO 


<r 


<r 


in 


00 


<f 


ON 


ro 


un 


ON 


un 


ON 


un 


X 


in 


— 1 


m 


o 


o 


CO 


cO 


CO 


o 


w • 4 


ON 


Mf 


wo 


ON 


in 


ON 


m 


+ 


oo 




00 


ON 


O 


—4 


i—4 


—4 


ON 


00 


NO 


00 




NO 




nO 




. — ( 






















• 


• 


• 


• 


• 


♦ 


• 


v-^ 










i— 4 


f— *4 




»' *' 1 
























































— 

CL 


r^ 


ON 






CM 


r>. 


o 


CM 


MD 


NO 


MD 


MD 


Mf 




O 


CO 


MD 


B 


i£> 


Mf 


<r 




CM 


un 


CO 


vO 


MD 


Mf 


CO 




un 


NO 


un 


CM 

ry*\ 


CO 


> ' 


CM 


vO 


CM 


o 


CM 


00 


00 


MO 


ON 


NO 


On 


O 


O 


i—< 


ON 


vJJ 


00 


* 
































# 


• 


P 


’~ m> 




— ( 


i—l 


r—i 














—4 


*—4 












O 


m 


r*. 


CO 


<T 


—4 


m 


ON 


ON 


<f 


CO 


00 


<r 


00 


NO 


co 


nO 


CV 


ON 


ro 


00 


—4 


00 




O 


ON 


CO 


ON 


CO 


« ■ 4 


— i 


ON 


co 


nO 


m0 


M. 


ON 


ON 


ON 


o 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


ON 


00 


ON 


ON 


On 










1—1 




























CM 






































o 


—1 


ON 


nO 


ON 


i—4 


—4 


MD 


MS' 


CM 


MD 


co 


00 


o 


O 


o 




X 


CO 


MO 


MO 


c- 


ON 


un 


CM 




o 


CO 


c*. 


r—A 


un 




O’ 




o 


+ 


ON 


ON 


ON 


ON 


ON 


o 


—4 


— I 


—4 


CM 


— 1 


f— 4 


— i 


i—4 


i—i 


p^4 


r— 1 


r—k 




































V f 




















—4 




i—4 


—4 


•—4 




r-A 


— *4 


<1> 


o 


o 


O 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


O 


6 


o 


ro 


o 


CO 


O 


CO 


o 


CO 


o 


CO 


o 


CO 


o 


co 


o 


CO 


o 


-.-4 


in 


un 


nO 


nO 


c* 


f"*. 


00 


oo 


ON 


ON 


o 


O 


c— 4 


»—4 


CM 


CM 


co 


•J 


— < 1 


— • 




— • 


1—4 




•—4 


— f 


*— 4 


—i 


CM 


CM 


CM 


CM 


CM 


CM 


CM 



•44 










' 







I 



